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Abstract 



D , We study by numerical simulation a lattice Yukawa model with naive fermions at inter- 

^ ' mediate values of the Yukawa coupling constant y when the nearest neighbour coupling k 
of the scalar field $ is very weakly ferromagnetic (k ~ 0) or even antiferromagnetic {k < 0) 
^ ■ and the nonvanishing value of ($) is generated by the Yukawa interaction. The renormalized 
Yukawa coupling j/^ achieves here its maximal value and this ^/-region is thus of particular 
importance for lattice investigations of strong Yukawa interaction. However, here the scalar 
field propagators have a very complex structure caused by fermion loop corrections and by 
the proximity of phases with antiferromagnetic properties. We develop methods for analyz- 
ing these propagators and for extracting the physical observables. We find that going into 
the negative k region, the scalar field renormalization constant becomes small and does 
not seem to exceed the unitarity bound, making the existence of a nontrivial fixed point in 
the investigated Yukawa model quite unlikely. 
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1 Introduction 



The scalar field self-coupling and the Yukawa couplings in the electroweak theory are be- 
lieved to vanish in the limit of infinite cut-off, as suggested by the signs of the perturbative 
/9-functions. But the confirmation of this "triviality" of these couplings and a reliable deter- 
mination of its consequences require nonperturbative methods, because one has to control a 
very difficult regime - when the renormalized couplings have maximal possible values which 
could, in principle, be large. A use of nonperturbative lattice methods for these purposes is 
desirable. 

In the pure scalar sector of electroweak interactions these methods have been successful 
in estimating the upper bound for the Higgs boson mass in the approximation neglecting the 
gauge and fermion fields [ip^li]- They also indicate that the renormalized quartic coupling 
is lower than the upper bound obtained from the tree level unitarity. This result makes a 
strongly interacting Higgs sector rather improbable and explains the quantitative agreement 
with calculations using the renormalized perturbation theory. Weakly coupled gauge fields 
do not seem to have any unexpected effects on the theory and can thus be treated 
perturbatively. 

Possible effects of a strong Yukawa coupling remain relatively unexplored, however. The 
most important phenomenological questions are: 

• What would be the effect of a strong Yukawa coupling on the upper bound for the 
Higgs boson mass 0? 

• What is the lower bound on the Higgs boson mass which follows from the vacuum 
stability requirement in the case of a strong Yukawa coupling H? 

• How strong can the Yukawa coupling actually be? What is the maximal fermion mass 
which can be generated by a Yukawa interaction |P? 

From the general point of view of quantum field theory the investigations of models with 
strong Yukawa coupling attempt to elucidate the following issues: 

• Are the Yukawa theories in 4 dimensions really trivial or do some nontrivial fixed points 
exist? 

• If they are trivial, how far and in what form Yukawa theories can be used as effective 
quantum field theories? 

• What are their relations to the purely fermionic theories with four-fermion coupling 

Recently numerous explorations of various lattice Yukawa models have been performed 
(for a recent review see ^^). In the following discussion, only the models with the 



so-called lattice parametrization and single-site Yukawa coupling are considered. The studies 
reveal the existence of three qualitatively different regions of the bare Yukawa coupling y. 
If y is sufficiently small, the lattice Yukawa models behave according to the perturbation 
theory based on the quasi-classical picture of the spontaneous symmetry breaking (SSB) 
in the scalar sector. For large y the fermions decouple in the continuum limit so that the 
physically interesting range of the values of y is restricted from above. 

In the intermediate region, the phase with nonvanishing scalar field expectation value 
($) exists even if the nearest neighbour coupling k of the scalar field is weak or even antifer- 
romagnetic (k < 0). The corresponding SSB is then generated by the Yukawa coupling. In 



particular, at k = the Yukawa models on the lattice correspond to the pure fermion theo- 
ries with multifermion couplings of the Nambu-Jona-Lasinio type [10|-|n| and one realizes 



that the SSB in such theories can be understood on the lattice as a special case of the SSB 
caused by the Yukawa interaction. 



In the light of the above-mentioned goals a very important observation |^ is that it 
is this intermediate region where the renormalized Yukawa coupling achieves its maximal 
values. Thus there is ample motivation to investigate SSB in Yukawa theories on the lattice 
in the regime where the Yukawa interaction is its driving mechanism and the standard quasi- 
classical approximation based on the scalar field potential is not adequate. 

Our recent numerical investigations of the SU(2)i®SU(2)K lattice Yukawa model with 



naive lattice fermions have shown ^ that in the region of intermediate values of y 
and negative k the scalar field propagators have a very complex form. As it turns out, it 
has two reasons: Firstly, the Yukawa coupling causes here sizeable fermion loop corrections 
to the scalar propagator. These effects can actually be used to estimate the value of the 
renormalized Yukawa coupling. Secondly, the antiferromagnetic ordering tendency of the 
negative scalar field coupling k competes with the ferromagnetic ordering effect of the Yukawa 
interaction. This shows up at large momenta and is therefore a lattice artefact which, 
however, has to be taken into account particularly on small lattices in the analysis of the 
scalar propagators. 

In this paper we develop reliable methods of analysis of numerical data in Yukawa models 
and present some results for the renormalized Yukawa and scalar quartic couplings. The main 
improvements over previous investigations of a similar type are threefold: 

(i) We have been able to analyze the scalar propagators in a very satisfactory way by 
including the 1-fermion loop contribution and as a result the Goldstone wave function 
renormalization constant is now determined quite precisely. 

(ii) To gain experience with finite-size effects we have varied the lattice size upto 12^16 
and have found the main effect to be the Goldstone boson dictated dependence. 

(iii) We have also covered nearly the full range of bare parameters which includes the 
maximum possible bare Yukawa coupling, staying within the lattice artefact-free region 
of the coupling parameter space. 

However, the fermion number in our simple model with naive fermions is too large. This 
probably causes the most serious problem we have found, namely that the fermion mass stays 
much below the scalar mass even for the largest renormalized Yukawa coupling. Therefore 
phenomenologically relevant large scale investigations of Yukawa theories should be per- 
formed in more realistic models. 

After defining the model in sect. 2 we describe in sect. 3 the spectrum in various phases 
and point out the appearance of the "staggered" scalar states to be interpreted as lattice 
artefacts occuring on a hypercubic lattice in the vicinity of the phases with antiferromagnetic 
ordering. In sect. 4 we discuss the physical meaning of the SSB occuring at negative k. 
The complex structure of the scalar propagators and the method of their analysis on finite 
lattices are described in sects. 5-7. In sect. 8 the fermion and scalar masses are discussed. 
Some results for the renormalized Yukawa coupling and the renormalized scalar field quartic 
coupling are presented in sects. 9 and 10, respectively. We summarize our results and 
conclude in sect. 11. Also a brief appendix elucidating the properties of models with both 
the ferromagnetic and antiferromagnetic orderings is included. 



2 The SU(2)/^(8)SU(2)/? model with naive fermions 



Our starting point is a fermion-scalar system in the continuum, which is defined in euchdean 
space-time by the action 
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(2.1) 



Here $o is a 4-component scalar field with a quartic self-coupling and \l/o is a fermionic SU(2) 
doublet field coupled to the scalar field by an Yukawa interaction with the bare Yukawa 
coupling parameter y^. and are the right- and left-handed chiral projectors. Because 
both fermions in the doublet couple with the same strength j/o to the scalar field, the action 
has a global SU(2)r fiavour symmetry. Together with the global SU(2)£, symmetry, which 
would turn into a local symmetry when gauge fields are included, the action is invariant 
under the global chiral 811(2)^(8)811(2)^^ transformations 
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(2.3) 
(2.4) 



where VtL,R e '&^{'^)l,r- 

We regularize the model by introducing a 4-dimensional hypercubic lattice with the lattice 
spacing a. A simple possibility is to keep the continuum parametrization and just to replace 
the derivatives by the lattice differences in the action (|2.1| ). But for a study of the largest 
renormalized Yukawa coupling it turns out to be very important to rescale the fields 



and reparametrize the coupling parameters 

(amo) = 8 , 
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(2.6) 



thus ending up with the model in the lattice parametrization defined by the action 
S = -^Y: ^Tr {$t$^^. + $t^^$^.} + ^ iTr j$t$, + A ($t$, - 



X/l 



We note that all the fields and the parameters in this expression are dimensionless. In the 
following we set the bare quartic coupling A of the scalar field to infinity, which leads to a 
freezing of the radial mode of the scalar field, |Tr = 1, so that can be represented by 
an 8U(2) matrix. The hopping parameter of the scalar field k and the bare Yukawa coupling 
y are left as free parameters. As usual in statistical mechanics, in the action ( p.7|) also 
the values k < are admissible, though the relationships ( p.5|) and ( |2.6|) to the continuum 
parametrization are not defined for negative k. 

For the sake of simplicity we are using the naive lattice fermions and the model actually 
involves 16 degenerate fermion doublets. There are some approaches trying to achieve a 



more realistic fermion spectrum preserving the chiral symmetry of the action. So far none of 



these approaches have proved satisfactory (see e.g. and for recent reviews |^ ^). An 
analytic treatment of the naive model by the 1/N expansion is not feasible in the investigated 
limit A = 00. In order to use the Hybrid Monte Carlo algorithm in the numerical simulations 
we further have to double the number of fermions by squaring the fermion determinant. Thus 
we are actually simulating the model with 32 fermion doublets expecting that qualitative 
aspects of Yukawa models are not changed by a large fermion number. 

One of such properties, found in many lattice Yukawa models with different symmetries 



and numbers of fermions |jT6|-||2^, is the occurence of SSB even at negative values of the 
hopping parameter k. Here we face the somewhat puzzling fact that the relation ( p.6|) 
between the coupling parameters of the continuum (|2.1|) and lattice parametrizations ( p.7| ) 
breaks down for negative k. As long as physical (renormalized) quantities are considered, 
the region /t < can have a well defined physical meaning. The physically interesting region 
of the parameter space seems to be extended in the lattice parametrization and the use of 
this parametrization is thus crucial. 

This observation motivates definitions of renormalized quantities which do not use the 
field $0 nor the transformations ( [^.5D and ( |2.6| ) and are thus applicable also for k < 0: The 



scalar field in lattice parametrization is renormalized in the phase with SSB as 

= ^ , (2.8) 



where Z^^ is the wave function renormalization constant of the Goldstone components of 
the $-field propagator. For k > this is equivalent to the renormalized scalar field in the 
continuum parametrization 

= ^ (2.9) 



when the dimensions of the fields are taken into account, 

= a <l>o,fl(x) , (2.10) 

and 

Zo,^ = 2kZ^. (2.11) 

The vacuum expectation value of the renormalized scalar field in lattice units is obtained 
from the magnetization ($) 

avR = ^ (2.12) 

and Vji is the vacuum expectation value in physical units. The renormalized couplings in the 
broken phase may be defined as 

yR- — -^^z^, (2.13) 

and 

where mp and m„ are the fermion mass and the a boson mass, respectively, in physical 
units. 



3 Spectrum and continuum limit at negative k 



The phase diagram of the 811(2)^(8)811(2)/? model (|2.7| ) with naive fermions shown in fig. 1 
[ p!9| , ^ consists of several phases and phase regions with ferromagnetic (FM), paramagnetic 
(PM) , antiferromagnetic (AM) and ferrimagnetic (FI) ordering of the scalar field. In addition 
to the weak coupling phases FM(W), PMW and AM(W)[] we find at large values of y the 
strong coupling phases FM(8), PM8 and AM(8) with nonperturbative behaviour of several 
fermionic observables. A similar phase structure was also observed in fermion-scalar models 
with different symmetry groups and/or different formulations of fermions on the lattice ||T6| - 
2^ provided the lattice parametrization and single-site Yukawa coupling is used. Common 



features of all these models are (i) the existence of weak and strong Yukawa coupling phases 
and phase regions, (ii) the continuation, at intermediate values of y, of the FM phase into the 
negative k region and (iii) the existence of points where several phase transition lines meet. 
The width of the funnel in the phase diagram filled by the FM and FI phases becomes smaller 
if the number of fermions is decreased (compare, e.g., the phase diagrams in refs. Q) 
in accordance with mean field calculations ||21[| . 

The fermion mass obeys within the error bars the tree level relation amp = y ($) 



in the FM(W) phase and is zero in the PMW phase where ($) = 0. It increases when the 
FM(8)-PM8 phase transition is approached from the FM(8) side, does not feel this phase 
transition and continues to increase when k is lowered within the PM8 phase. 

The scalar spectrum is quite different in different phases. In the four phases around the 
point A, where the PMW, AM(W), FM(W) and FI phases meet, the spectrum is indicated 
in fig. 2. In the phases PMW and AM(W) with zero magnetization ($) there exists a scalar 
$ particle quadruplet of degenerate mass am$. On the other hand, in the phases FM(W) 
and FI with nonzero magnetization ($) there exist three Goldstone bosons called vr and one 
massive a boson with masses am7r = and arricr, respectively. 

In the K < region the scalar propagators show near the momentum (tt, vr, vr, vr) the 
presence of further scalar states which we call staggered states in the following. This effect 
shows up feebly already at small positive k and as k is lowered it results in a gradually 
increasing curvature in the scalar propagators at large momenta. The existence of such 



states is obvious at y = 0, as here the symmetry of the action (|2.7] ) under the transformation 

^ $f = (_l)^'l+^2+X3+X4^^^ K-^-K (3.1) 

implies the presence of a quadruplet in the PM phase as well as of three vr''* and 
one cr''* states in the AM phase with nonvanishing staggered magnetization 
^$■5*^ = (^^(— l)^i+^2+a;3+a;4<|,^^^ They are visible at low momenta in the 2-point func- 
tion of the field "l*^*. We denote the corresponding masses am|*, am^ and am^J, respectively. 
For < y < oo the symmetry ( p.!] ) is broken explicitly, and these states can appear on the 
lattice simultaneously with the usual $ or tt and a bosons, as indicated in fig. 2. 

The question is, which particles will remain in the various possible continuum limits, or 
in the large cut-off limits if the theory is an effective one (several possibilities are discussed 
in ref. |]2^). We want of course to recover in the continuum limit the fermions and the 
usual bosons simultaneously. This includes the renormalized vacuum expectation value of 
the scalar field avR = which is proportional to the gauge boson mass amw This is 

possible only in the FM(W) phase in the scaling region of the FM(W)-PMW phase transition. 
Here the model is of physical relevance from the point of view of the electroweak theory. 



""^For simplicity we use the nomenclature "phases" also to denote phase regions with weak (W) or strong 
(S) bare Yukawa coupling and include a bracket in the abbreviation. 



As am|* and am^* vanish only on the critical lines where is vanishing, the staggered 
states could remain in the continuum limit simultaneously with the fermions and the usual 
scalars only if the limit is taken at the point A of the phase diagram in fig. 1. The scaling 
region of this point should therefore for physical reasons probably be avoided. The staggered 
states are lattice artefacts because they depend substantially on the lattice geometry. But 
they can have drastic effects on finite lattices and thus have to be taken into consideration 
in the analysis of the numerical data for the scalar propagator. 

4 SSB generated by Yukawa coupling and the relation 
to four-fermion coupling 

The negative k region has no analogue in the continuum parametrization because the trans- 
formation equations ( |2.6| ) are not defined for k < 0. At first glance it seems to be awkward 
and the question immediately arises whether a sensible continuum limit can be obtained in 
this region. The fact that for small positive k and for negative k there is a broken symmetry 
phase FM(W) at large enough Yukawa coupling suggests that the Yukawa coupling must 



be the driving force for SSB in that region (see also ref. |2^). This is obviously outside 
the regime of usual perturbation theory in the continuum. It is therefore very important to 
investigate this region if one is looking for nonperturbative effects of the Yukawa coupling. 

Of course, the AM and FI phases appearing only for k < are probably lattice artefacts 
as they depend very much on the lattice geometry. But the scaling regions of the FM(W) 
and PMW phases down to the point A (fig. 1) are worth of consideration. 

The Osterwalder-Schrader (OS) reflection positivity can be proven presumably only for 
K > 0. However, it is a sufficient, but not a necessary condition for unitarity, so that unitarity 
can still hold. In the numerical computations of the propagators of the theory, our failure 
to detect any state with a negative norm is assuring. Furthermore, as demonstrated later in 
this article and also in ref. , the measured values for all the masses and the renormalized 



couplings continue analytically from positive to negative k across k = 0. 

The problems with the transition from k > to k < seem to exist only on the level of 
bare parameters in the continuum parametrization. One should consider only the renormal- 
ized quantities. If one knew the renormalized running coupling y(/i) at all momentum scales 
H one could define also a sensible bare Yukawa coupling |/b as 

Z/B = ~ • (4.1) 

For small values of Uq and g^, usual perturbation theory is valid and ub would not differ 
very much from the coupling parameter uq. At larger values of and in particular in 
the negative k region the parameter does not any longer have to be close to ys- This so- 
defined bare coupling ys would also have in general no simple relation to the bare parameters 
ft and y. Fig. 3 schematically illustrates how yB and y^ may split up from each other as 
the nonperturbative region is entered along some line of constant physics specified, e.g., by 
— f^phy) — const, where /^phy is a physical scale. 
The stage is therefore well set to plunge into the negative k region with three issues in 
mind: 



(i) Is there a nonperturbative fixed point? 

(ii) Even if the theory is trivial, can the couplings be strong at a reasonably large cut-off? 



Our previous investigation |^ with more naive methods of analysis of the numerical data 
has already produced tentatively negative answers to these two questions which we want to 
confirm in this paper. 

(iii) In any case, it is necessary to find out how far the observables in the k < region 
differ quantitatively from the k > region, for determination of bounds on renormalized 
couplings. 

In the rest of this section we want to point out the connection of a Yukawa theory with 
a Nambu-Jona-Lasinio (NJL) type model. Integrating out the scalar fields in the partition 



function defined by the action ( p.7|) with k = X = produces obviously a pure fermionic 
theory with local (on-site) four-fermion coupling of strength ^y"^ - the NJL model on the 
lattice. The scalar field is equivalent to the "auxiliary field" used in the context of the 
four-fermion coupling [0, |lT], [3^]. Also for A > the effective fermion interaction is local 
at K = and corresponds thus to a multifermion interaction discussed, e.g., in ref. [|12|. 
Thus we conclude that theories of the NJL type are special cases (k = 0) of the Yukawa 
models in the lattice parametrization. But in terms of the renormalized theory, as we have 
discussed above, k = is not singular and the same qualitative physics is obtained for the 
whole FM(W)-PMW scaling region down to the point A unless there is a nonperturbative 



fixed point somewhere. Recent work |]T3|, using expansion also shows an equivalence 
between Yukawa models and generalized NJL models. 

Thus we achieve a unification of concepts and language of the SSB: the Yukawa models 
treated nonperturbatively using the lattice formulation ( p.7|) embrace both 

• the classical Higgs mechanism, in which the SSB is understood in terms of the quasi- 
classical approximation for the effective potential of the scalar field and perturbation 
theory (e.g. the region ?/ -C 1 in fig. 1 for any A), and 

• the NJL type mechanism, operating at small (and possibly negative) k, for relatively 
large values of y, which has to be treated by nonperturbative techniques such as the 

expansion. 

Yukawa theories provide a gradual transition between these mechanisms. It seems, therefore, 
possible to formulate the SSB in the standard model in terms of the NJL mechanism |]T0|, |Tl|, 
[T^ . However, in the light of the above discussion, the distinctions between an elementary 
and a composite scalar, an auxiliary and a dynamical scalar field and between dynamical 
symmetry breaking and the usual Higgs mechanism do not seem important. 



5 Properties of the scalar propagators 

In our numerical simulation we consider V = L^T lattices with periodic boundary conditions 
for the scalar fields. Fermionic fields are periodic in space and antiperiodic in euclidean time 
directions. 

In the symmetric (PM) phase the bosonic spectrum contains the $ quadruplet of 
mass am$. The corresponding propagator in the momentum space is 

G^iap) = |Tr{$t$Jexp(zap(x-y))^ . (5.1) 



Neglecting the instability of $ at the present precision level, the renormalized mass am$ and 
the wave function renormalization constant Z$ can be defined by means of the limit 



G<i.(ap)|p2^o = 7 , -^2 > (5-2) 



where the quantity ap = 2J2^i^ — cos(ap^)) is the dimensionless lattice equivalent of the 
momentum square in the continuum. 

In the broken (FM) phase it is useful to introduce the following notation for the scalar 
field 

3 

$ = al + z^TT^' r^' . (5.3) 



Here t\ j = 1, . . . , 3 are the usual Pauli matrices and ycr^ + J2j=i{'^xy = 1. The compo- 
nents are chosen such that the magnetization is given by ($) = (a). Then the longitudinal 
component is associated with the massive a boson whereas the transverse ones with the three 
massless Goldstone bosons vr. On the lattice the vr and a propagators in the momentum space 
are defined by 

Gniap) = ^^5]II^x7r^exp(mp(x -y))^ , (5.4) 
Ga{ap) = /-^^a^aj^exp(zap(x - ?/))\ . (5.5) 

The only asymptotic states in the FM phase are the massless vr bosons. Therefore the wave 
function renormalization constant for the scalar field is defined through the following 
limit of the vr propagator, 

G^{ap)\p2^Q = . (5.6) 
ap 



Using the so-defined Z^^ the renormalized field expectation value vr is then given by eq. (|2.12| ). 



Again at our present precision level it is presumably sufficient to define the renormalized mass 
arria of the unstable a particle by the relation 

G^(ap)|p2_,o = 7 ^2 • (5-7) 

[am^)^ + ap 

Another point to note is that in the pure $^ theory the renormalized mass defined this way 
is very close to the physical mass 0, |^. 

In a finite system no spontaneous breakdown of the symmetry can occur and during a 
Monte Carlo simulation in the broken phase the system drifts through the set of degenerate 
ground states. This causes a vanishing of noninvariant observables like ($). To compensate 
for this drift, each scalar field configuration is rotated so that ^ J2x = ^ J2x '^x- In the 
pure theory the rotation technique provides a very good approximation of the infinite 
volume values of the noninvariant quantities 0. 

A lot is known about the properties of the scalar propagators in the pure $^ theory 
which is the limiting case of the model ( p.7|) both for ?/ = and y = +oo (in the latter 
case fermions become infinitely heavy and decouple completely from the particle spectrum) . 
When plotting the inverse propagators G5^(ap), G~^{ap) and G~^{ap) as functions of ap^ one 
finds for all propagators and for all possible values of ap^ a straight line behaviour confirming 
the analysis of the data in terms of free scalar propagators. From straight line fits to the 



inverse propagator data, using the relations (|5.2|) , ( f).Qf i and (|5.71 ), one can determine the 
wave function renormahzation constants and the renormahzed masses on a finite lattice. 
In the Yukawa model we have determined for various values of k and y the momentum 



space propagator (|5.1| ) in the PMW and PMS phases and the propagators ( p.4[ ) and ( p. 51 ) 



in the FM phase close to the critical lines FM(W)-PMW and FM(S)-PMS (see fig. 1). For 
very small and very large values of the Yukawa coupling y the results for the propagators are 
very close to those obtained in the $^ theory, i. e., when plotting the inverse propagators as 
functions of ap^ we find approximately a straight line behaviour. However, when entering the 
intermediate Yukawa coupling region and lowering k, where fermions have a strong feedback 
on the scalar sector and the staggered scalar states (sect. 3) become visible, deviations from 
the free propagator behaviour at larger values of api^ are observed. These deviations become 
more and more pronounced when approaching the multicritical points A and B. 

In fig. 4 we display for several points in the FM phase the inverse Goldstone propagator 
G~^{ap) as a function of the quantity ap^. The figures in the left column were obtained 
at 3 points in the vicinity of the FM(W)-PMW phase transition whereas the figures in the 
right column correspond to 3 points in the vicinity of the FM(S)-PMS phase transition. The 
lowest figures were obtained very close to the multicritical points A and B respectively. The 
figures show that there are three kinds of deviations from a free propagator behaviour: 

1. The formation of a second pole in G.,t{ap) in the corner of the Brillouin zone with 
ap = {it, it, tt, vr) when approaching the points A or B within the FM phase. This effect 
is caused by the staggered states which are present on the lattice also in the FM 
phase (see sect. 3). 

2. The appearance of dips in around the momenta ap^ = 4, 8, 12 (corresponding to 
apf^ = or tt) in the weak coupling region. These dips are already visible at positive k 
as can be seen from the first figure in the left column. 

3. At small ap the inverse propagator has in the weak Yukawa coupling region a curvature, 
which plays a significant role in the data analysis. 

Similar structures were also discovered for the propagator Ga{ap) in the FM phase and the 
propagator G^{ap) in the PMW and PMS phases near the FM(W)-PMW and FM(S)-PMS 
phase transitions. As we shall now discuss, the first two effects are actually lattice artefacts, 
dependent on the geometry of the lattice, but they have to be understood quantitatively 
in order to extract the physically relevant quantities from the scalar propagator data. The 
third effect is physical. 



6 Staggered scalar states in the FM and PM phases 

For an understanding of the two pole phenomenon it is useful to discuss first the situation 
in the pure theory which is found in the limiting cases y = and y = +oo. We use the a 
propagator G„{ap, k) as an example and indicate for a moment explicitly the ^-dependence 
of the propagator. Using the transformation ( p.l|) we find the relation: 

G^{ap, —k) = G^{ap + ail, k) (6.1) 

where vr = (tt, tt, tt, tt). As the propagator G„{ap,K,) in the pure scalar theory can for all 
momenta ap be well described by the free scalar Ansatz eq. (|5.7] ), G^^ap, —k) has for k > Hc 



the form 

GAap, -n) = ^ , = , ^2 (6-2) 

{am,Y + a{p + tt) [am^^? + 16 - ap 

with m^* = m„ and = Z„. For am^ = the propagator G(j{ap, —k) has thus a pole at the 
momentum p = fc. Similar relations can be obtained also for the propagators and Gt^. 
Obviously, at no k the normal and staggered states appear simultaneously as their masses 
are small only close to the FM-PM and PM-AM phase transitions, respectively, which are 
distant. 

For < y < oo the symmetry does not hold any more. Nevertheless, for small and 
large values of y the spectrum is very similar to the pure "l"^ theory. But when the phase 
transition lines approach each other and finally meet at the points A and B we expect that 
both the normal and the staggered state masses are small simultaneously. In particular, 
am I* can be small also in the FM phase. In the FM(S) phase it is therefore reasonable to 
try to fit the numerical results for G^{ap) by the two pole Ansatz 

^'^M = 7 + . ^2 • (6-3) 

[amcr) + ap (am|, + lb — ap 

Analogous expressions have been applied also for the propagators G-^^ap) and G<!,{ap) in the 
FM and PM phases respectively. In fig. 5 we show as an example a fit to the Goldstone 



propagator in the FM(S) phase, where the fit Ansatz for G^(ap) is given by eq. (3^) with 
arrij, = and Z^. replaced by Z^. Fig. 5 shows that the data are described by the two pole 
Ansatz very well. We furthermore expect that the staggered mass am|f and the wave function 
renormalization constant Z^ obtained from the fits to G^^lak) and Gfj{ap) should agree. 
This expectation is indeed confirmed by the numerical results, for example at k = —0.65 
and y = 1.8 the values are am|* = 1.47(7) / 1.37(4) and Z|* = 1.41(7) / 1.45(4) from 
the (t/tt propagators. As expected, the mass am|* does approach zero when the line AB is 
approached. It should be also stressed that both terms in the Ansatz (|6^ ) have positive 
sign, so that the second term has the usual form of a pole after shifting the momentum by 
(vr, TT, TT, tt) and cannot be interpreted clS db ghost. 

The same formulae describe also the two pole structure of the scalar propagator in the 
vicinity of the point A. Here, however, a more elaborated Ansatz has to be developed in 
order to take into account simultaneously also the overlaid finer structure caused by the 
fermion loop corrections. It is the subject of the next section. 

7 One fermion loop contribution to the scalar propa- 
gators 

The other two features making the scalar propagators different from a free propagator are 
the appearance of dips at momenta ap^ = 4, 8, 12 and the curvature at small ap^. They occur 
in the weak coupling phase regions FM(W) and PMW where the fermion masses scale. 

According to the definitions of the quantities Z^, Z^- and am^^ in eqs. ( |5.6| ) and ( p. 71 ) the 
scalar propagators have to be analyzed in the limit ap'^ — > 0. However, on small lattices with 
periodic boundary conditions the smallest nonvanishing momentum is ap = (Jot T > L) 
which is quite large, e.g., ap = 0.79 for T = 8. In the pure theory this is not a serious 
problem as the inverse scalar propagator can be fitted with a straight line up to ap^ = 0(10) 



In our Yukawa model the situation is much less favourable: in addition to the dips 
occuring at a/p^ = 4, 8 and 12 - for which one might argue that they should simply be 
ignored as the analysis has to be restricted to the smallest momenta - we have to face 
the more serious problem of a significant curvature at small ap^. The following further 
observation makes the situation look even worse: Increasing the T-extent of the lattice - 
which is the cheapest possibility to have small momenta - e.g. from T = 6 up to T = 46, 
we still have not found an onset of a linear a/p^ dependence. We conclude from this that for 
sizeable Yukawa coupling an application of the free particle Ansatz for the scalar propagators 
on finite lattices produces uncontrollable systematic errors. 

Therefore, we have developed a more sophisticated fit Ansatz for the scalar propagator 
based on the 1-fermion-loop contribution to the self-energy of the vr or a bosons. The 
justification for including only the fermion loop comes from the experience in the pure scalar 
sector of the model where the scalar loop contributions do not change the linear op^-shape 
of the propagators but only lead to wave function and mass renormalizations. On the 
other hand, as will be shown below, in the case of naive lattice fermions the 1-fermion-loop 
contribution will cause deviations from the linear a/p -dependence just of the form observed 
in the data. 

Let us discuss the example of the Goldstone boson propagator. On finite lattices we may 
write 

G~\{ap) = [{am^^Lf + ap^ - T.^,L{ap; amp^L)] , (7.1) 

where denotes the 1-fermion-loop contribution to the self-energy of the Goldstone boson 
in the renormalized perturbation theory. For the purpose of this section the subscript L 
points out the possible dependence of various quantities on the spatial lattice size L. For 
instance with regard to the Goldstone bosons we take the finite-size of the lattice into account 
in a naive but simple way - allowing for a finite mass of the Goldstone boson airiT^^L (see 
e.g. also [0). (A treatment based on chiral perturbation theory, successful in the pure 
theory 0, is in the complex situation with light fermions presumably not applicable.) In 
the FM phase where the fermions are massive we impose the two necessary normalization 
conditions on Sj^ in the infinite volume at momentum zero, 

Svr,oo(ap; ami7'^oo)|p=o = 0' 

= 0. (7.2) 



(ap; amp 



dap 



p=0 



With this normalization approaches in the thermodynamic limit the form used for the 
definition of Z^^, eq. (p.6|), provided (am^r.L)^ — as L ^ cx3. 

The not-yet-normalized ^ calculated from the corresponding Feynman diagram on the 
lattice is 

4 f 1 ^ 1 1 

E'^iap; amF,L) = (-l)T^I]Tr <^ (^752/fl) tttt^ (^52//?) 



i^{k) + amp p i^{k — ap) + amp p J 

(7.3) 

This equation defines the renormalized Yukawa coupling in terms of the 3-point function. 
Of course, could in principle differ from i/fi defined in eq. ( |2.13| ). In the above s^{k) = 
sin kfj^, the factor (—1) comes from the fermion loop, the trace is to be taken over the 
Dirac indices. The SU(2) trace together with another factor of 2 from the HMC doubling 
(see sect. 2) results in the factor 4 whereas the standard fermion doubling is taken into 
account automatically. The sum over the loop momentum k runs over the set of momenta 



corresponding to a finite Li^T lattice with periodic boundary conditions in the space direction 
and antiperiodic boundary conditions in the time direction, 




rij = 0, 1,2, ...L- 1 
^4 = 0, 1,2,...T - 1 . 



J = 1,2,3 , 



(7.4) 



After some simplifications we get 



Tl^ L^ap; amF,L) 



16 ^ {amF,LY + s{k)s{k — ap) 



'""L^T ^ [(am^,z.)2 + s\l^)\ [(am^,^.)^ + s\k - ap)] 



(7.5) 



where we introduce the notation ^ for the (unnormalized) lattice integral after factorizing 
out y\. In order to recognize the 1-fermion-loop as the reason for the deviations in the scalar 
propagators it is instructive to discuss the dependence of this expression on the external 
momentum ap. For ap = (0, 0, 0, 0) the denominator has minima when the loop momentum 
components are fc^ = or tt, corresponding to processes involving the physical fermion 
and its antifermion or some doubler fermion with its own antidoubler. Thus a considerable 
contribution of to G-,^ for small momenta ap can be expected, causing the observed strong 
curvature of . In addition, also peaks when some components ap^ = tt and others 
are zero. Then the kinematics allows such intermediate states to be excited which involve 
e.g. the physical fermion and the antidoubler of momentum ap, or any other appropriate 
pair of doubler and antidoubler whose respective positions of poles differ just by ap. This is 
precisely the reason for the dips seen in the inverse scalar propagators near ap^ = 4, 8, 12. 

To fit the propagator data using the Ansatz ( |7.1| ) we have to perform the following steps: 
First, T^T^^i^ is normalized as 



so that it satisfies the conditions (|7.2| ) in the infinite volume limit. Here the fermion mass 
mp^L on the given finite lattice is taken from the standard fit to the fermion propagator data. 
But note that this normalization also requires an estimate of the fermion mass in infinite 
volume amp^oo (various attempts to normalize S,r,L using only finite volume quantities did 
not work). In the spontaneously broken phase FM, the major part of the finite-size effects is 
expected to be due to the massless Goldstone bosons leading to a volume dependence linear 
in 1/L^. So, checking this dependence and then extrapolating amp^i to anip^oo requires, at 
a given (k, ?/)-point, simulations on a sequence of at least three lattices. We have performed 
runs on lattices L^16 with L = 6, 8, 10 and 12 with the result that as long as amp^i itself is 
not too small the agreement with a linear 1/LP' dependence indeed allows an extrapolation 
to amp^oQ (see next section). 

Analogous to eq. ( [7.5|) we define the (normalized) lattice integral It,^l with yj^ factorized 
out. 



S^^i(ap; amp^L) 



E;^^(ap; amp^L) - K,ooi.(^P 



0; arrip^^) 




(7.6) 




(7.7) 



The TT propagator fit Ansatz is then 




with the free parameters am-i^^L, and yn (we note that through the normahzation condi- 
tions ( [7.2|) this expression also depends on mp^^). 

Before describing the results in the next sections let us discuss the quality of the fit. The 
full Ansatz we use is the superposition of the propagator with the pole at ap = (0, 0, 0, 0) 
including the 1-fermion-loop contribution and the staggered propagator with the pole at 
ap = (tt, TT, TT, tt) as explained in sect. 6 above. This Ansatz is able to describe the scalar 
propagator data in the complete interval < op^ < 16. In particular, the curvature at small 
momenta and also every detail of the peculiar structures near ap^ = 4, 8, 12 are perfectly 
reproduced. This is demonstrated in figs. 6 and 7 for two typical examples of the fits, one at 
small positive k and one in the close vicinity of point A where the scalar propagators have 
the most complex form. 

It should be stressed again that the dips at momenta ap^ = 4, 8, 12 are caused by the 
presence of the doubler fermions and hence should at least look different, if not absent, in 
models without them. However, in all models the curvature at small momenta will appear 
for sufficiently strong Yukawa coupling and the second pole at op^ = 16 will be present near 
phases with antiferromagnetic ordering. 

8 Fermion and scalar masses at small positive k and 
negative k 

We have been able to determine reliably both amp and arrirj simultaneously only for n > 0. 
Here ama is always greater than amp at least by a factor 3-6. This presumably is a 
consequence of having a large number (32) of fermion doublets. Estimating that a similar 
mass-ratio holds also for k < we have mostly performed calculations at points with very 
small fermion mass, amp ^ 0.1 — 0.3, in order to have the cj-boson mass at least smaller 
than 1. 

The determination of amp in this range of values requires only moderate statistics and 
can be reliably performed at k < even in the vicinity of the point A. An important 
condition is, however, that the long size T of the lattice L^T is at least 16, otherwise amp 
is spuriously small and T-dependent when determined from the fit by means of the free 
fermion propagator. The method of analysis and many results have been presented already 



in ref. ||2^. Here we would like to point out the large spatial volume dependence of amp. For 
amp ^ 0.3 the value can decrease by 40% when L increases from 6 to 10. Nevertheless, for 
amp ~ 0.2 the decrease is linear with so that one can tentatively extrapolate to L = oo. 
The observable ($) is easily measurable everywhere and also has a linear dependence. 
The dependence and the extrapolation of both the quantities are shown in fig. 8 (we 
show avR instead of ($) because Z.„ is L-independent as discussed in sect. 9). 

The a boson mass has been determined at several points for /t > on lattices of various 
sizes. The finite-size effects are compatible with the expected dependence but the large 
error bars prevent a verification. Nevertheless, we have used the same volume dependence 
to extrapolate am^ to L = cxd. 

There are two technical reasons making a reliable determination of am^ at negative k very 
difficult. Firstly, the number of iterations for the fermion matrix inversion needed for the field 
update increases drastically with decreasing k. As the determination of the a propagator 
requires much higher statistics than of the fermionic propagator, an accumulation of good 
data in the negative k region, in particular in the vicinity of the point A is prohibitively 
expensive. Secondly, the maximum of the curve G~^{ap) occurs already at rather small 
momenta (see fig. 7) and only a few data points of the propagator contain the information 



about ania and Zg-, the rest being dominated by the staggered state. Therefore we have not 
succeeded to determine rehably amo- for any of our data points at negative k. 

The largest renormahzed couphngs are expected on the boundary of the scahng region, 
i.e. relatively far from the critical line. However, at present we do not know the position of 
this boundary, actually not even its proper definition (e.g. how small amp and ania- should 
be in order that the lattice model can be used as an effective continuum theory). We are 
thus not able to extract upper bounds on masses from our results for renormahzed couphngs. 



9 The renormalized Yukawa coupling 

The excellent agreement of the fits with the MC data for the vr propagator both for positive 
and negative k allows to determine Z^^ reliably. In comparison to the usual determination of 
by a naive free particle fit to the smallest momentum used e.g. in our earlier publication 
[ p^ the present method yields more precise results (the error bars are reduced by factors 
3 - 10). In particular, they are now stable when the T-size of the lattice is varied, in 
spite of the strong curvature near = which is most clearly seen on large-T lattices. 
No L-dependence of has been found. The actual values are slightly smaller than found 



in ref. |^| confirming the conjecture in that paper that a simple linear fit to the smallest 
momentum on small lattices gives overestimated values of Z.,^. 

Some values of obtained in the vicinity of the FM(W)-PMW critical line are plotted 
in fig. 9. decreases strongly as k decreases and appears to vanish at the multicritical 
point A. 

The knowledge of Z-j^ allows us to determine, from the available very good data for airiF 
and ($) for various L, the renormalized Yukawa coupling by means of the definition (|2.13| ). 
Both amp and ($) are strongly L-dependent, but their ratio turns out to be practically L- 
independent. We use the observed linear dependence to extrapolate amp and avR to 
infinite volume obtaining y^i in the L ^ oo limit (see fig. 8) In fig. 10 we plot these results 
against the fermionic correlation length = l/amp at L = oo. The reason for not using the 
smaller = 1/ at L = cxd for this purpose is the fact that we do not know its values for 

K<0. 

The dotted curve in fig. 10, given by 



47r2 



' lna/i| 



is obtained by choosing infinite bare Yukawa coupling in the 1-loop formula for the running 



Yukawa coupling and identifying the scale fx to he nip. As pointed out in |]23|, this curve 
described quite well the results for at small positive k, and negative k. The dotted error 
bars associated with this curve in the figure indicate the range of values and the error bars 



of those former results in [23 



The dramatic reduction of the error bars in fig. 10 is mostly due to the refined analysis 
in this paper. It is now apparent that y^ is bounded by the dotted curve, suggesting 
applicability of 1-loop perturbation theory even close to the point A. In addition all the 
yji results are clearly below the s-wave tree level unitarity bound for 32 fermion doublets. 
This bound is indicated by the horizontal dashed line in fig. 10. These results suggest the 
conclusion that in our model there is no strong Yukawa coupling regime at least for C,f > 5. 
We remark that initial indications from numerical investigations in the mirror fermion model 
2^ are, however, different. The weakness of yn in our model in the negative k region 
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Table 1: Comparison of yji from the tree level definition l \2.1^) to yji from the fit to the 
Goldstone propagator (\7.<^ at some typical {K,y)-points. 

comes about as follows: the unrenormalized ratio amp/{^) actually increases strongly as k 
decreases, but this is compensated by the equally strong decrease of (fig. 9). 

The fermion loop correction allows to determine the renormalized Yukawa coupling 
defined in terms of the vertex function and extracted from the Goldstone propagator data 
by the Ansatz (|7.8| ). In table |l| we compare yn and yn for some typical (k, ?/)-points where 
we were able to extrapolate amp to L = oo in order to satisfy the normalization conditions 
( [7l2|) . The good agreement found between y^ and yn indicates that the analysis of the scalar 
propagator data by taking the 1-fermion loop corrections into account is adequate and that 
both definitions of the renormalized Yukawa coupling quantitatively agree. This further 
supports the conclusion that in our model with naive fermions the Yukawa coupling is not 
strong. 

Some caution is due, however. It could be that our linear extrapolation of amp to L = oo 
underestimates ^j?. Furthermore, as we do not know am^ in the vicinity of the point A, we 
cannot exclude that is as large as and that our data points there are actually deep in 
the scaling region. 

10 Fermion influence on the scalar mass bound 

To study the influence of heavy fermions on the triviality bound for the scalar mass it 
would be again most interesting to perform this analysis in the vicinity of the multicritical 
point A and the lack of reliable results for arricr there is deplorable. However, provided the 
renormalized Yukawa coupling does not attain in the negative k region values significantly 
larger than at k positive, as suggested by our results, we expect that it is sufficient to 
investigate the influence of the Yukawa coupling on the scalar sector in the positive k region. 

These considerations have motivated our relatively high statistics study of the cr-propaga- 
tor at K ~ 0. We have fixed y at the value y = 0.6, for which the critical k is Kc = 0.020(5). 
On lattices of three different spatial sizes, L'^16 with L = 6, 8, 10, we have accumulated 
about 4-5 thousand Hybrid Monte Carlo trajectories at three points k = 0.03, 0.04, 0.06. 



Both scalar propagators have been analyzed by means of the same Ansatz (|7.8| ) and the 
scalar mass arricr and have been determined. 

The results for the ratio ma/vR are displayed in fig. 11 as a function of the scalar correla- 
tion length S,a = l/aw^cr (open symbols). The finite-size dependence of am„ can be described 



to be linear in though the error bars on ama leave room for modification. The tentative 
extrapolation of the results to the infinite volume assuming a dependence is indicated 
by the full squares. This kind of plot is customary to extract a triviality upper bound for 
the scalar mass in the pure $^ theory. For comparison the data from this theory {y = 0) on 
a hypercubic lattice are also shown (full circles). 

The results in the Yukawa model on finite lattices approach the infinite volume results 
from below (different from the theory 0) and the extrapolated results are - within large 
error bars - consistent with the pure results. As ^o- ~ 1 - 3 and C,f is much larger than 
^o- for all points in fig. 11, one can expect that the edge of the scaling region is contained. 
So we observe no large influence of the Yukawa coupling on the Higgs mass upper bound in 
our model for k > 0. 



11 Summary and conclusions 

We have explored the region of the largest renormalized Yukawa coupling in a lattice Yukawa 
model with naive fermions in the broken symmetry phase. In this region the Yukawa in- 
teraction is the driving force of the spontaneous symmetry breaking overwhelming the very 
weak ferromagnetic or even antiferromagnetic nearest neighbour scalar field coupling. Such 
competing interactions, together with sizeable fermion loop corrections, result in a complex 
structure of the scalar propagators. We have demonstrated that this structure can be theo- 
retically understood and even utilized for an extraction of the renormalized Yukawa coupling 
yn from the scalar propagator data. The results for y^ showing that the Yukawa coupling 
is small in the limit of large cutoff are consistent with the triviality of this coupling. In the 
physically relevant FM(W) phase the lines of constant y^ seem to flow nearly parallel to 
the FM(W)-PMW phase transition and for y > 1 run out of the FM(W) phase, instead of 
flowing into some point on this phase transition. In particular the suspicious point A does 
not seem to be a nontrivial fixed point. The values of do not seem to exceed significantly 
those at K > 0, indicating that the k, < region of the FM(W) phase does not add much to 
the physical content of the model at k > 0. We expect these results to be generic for various 
lattice Yukawa models. 

Our results, looked at quantitatively, may however be specific to the chosen model with 
a large number of fermions (32 doublets). The renormalized Yukawa coupling stays below 
the tree level unitarity bound, thus being never strong. Correspondingly, no influence of 
the Yukawa interaction on the upper bound for the cr-mass could be detected. However, a 
word of caution is warranted because we have not localized reliably the edge of the scaling 
region, where the largest values of the renormalized couplings should actually be determined. 
Furthermore, the model suffers from a drawback caused by the large number of fermions: 
the fermion mass generated by the Yukawa interaction stays substantially lower than the a 
boson mass. This makes investigations in the scaling region difficult with two very different 
correlation lengths and is certainly not a generic feature of Yukawa models. There are now 
suggestions of Yukawa models with a small number (2 and 4) of fermion doublets . The 
methods developed in this article should be appropriate also in these models. Apart from 
quantitative differences, it would be interesting to see if any of the qualitative conclusions 
drawn in this paper change. 
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Appendix A: The metamagnetic Ising model 

Integrating out the fermionic fields leads to a scalar model with nonlocal interaction 
terms. In this appendix we want to discuss a simple type of spin model with a nearest- 
neighbour {nn)- and a next-to-nearest-neighbour (rm?T,)-interaction terms, which appears 
in the literature in the context of metamagnet^. In this model an overlap of the scaling 
regions associated with the normal and the staggered magnetizations occurs analogous to 
the vicinity of the points A and B in Yukawa models. 

The model is defined by the action 

g 

x,fi>0 x,fi>u 

where ax are Ising spins. 

In a mean field approximation the critical lines are given by p 



<n + <nn = ^^ising «nn > 0, > for the FM-PM transition line 

— -I- K^„„ = Kising f^nn < 0) '^nnn > ^^c AM-PM transition line, 

where the constant Kf^j^g is known from the 4-dimensional Ising model to be about 0.0748. 
At = the transition lines meet in a multicritical point. The phase diagram found in 
a numerical Monte Carlo simulation of the 4-dimensional model shown in fig. 12 has three 
phases: a ferromagnetic phase (FM) with (a) > 0, (o"'^*) = 0, a paramagnetic phase (PM) 
with (o") = 0, (c**) = and an antiferromagnetic phase (AM) with (a) = 0, (o"**) > 0. The 
phase transition lines separating the symmetric from the broken phases is of second order 
and the phase transition separating FM and AM is of first order. Except for the absence of 
the FI phase, this phase diagram is similar to the phase diagram of the Yukawa model at 
weak Yukawa coupling, the triple point being analogous to the point A in fig. 1. 

In the vicinity of the multicritical point we find by numerical simulation that the inverse 
propagator of the scalar field in momentum space has two poles, one of them at ap = 
(tt, tt, tt, tt) and some structure for intermediate momenta. An example in the ferromagnetic 
phase near the multicritical point is shown in fig. 13. This shape is well described by the 
inverse propagator of a gaussian model with the same (nn)- and (nnn) kinetic terms 

g 

Sp = -ann " cosA;^} - a„„„— ^{2 - cos(/c^ + K) - cos(/c^ - K)] + h. (A.2) 

The fit of the Monte-Carlo data with the function ( [A.2| ) plotted in fig. 13 demonstrates 
that one can describe the data very well in spite of a nonlinear dependence on ap^. If one 
includes further interaction terms, the structure of the propagator at intermediate values of 
aJp changes. 

^See footnote on p. 60 in |3^. 



Of course, integrating out the fermion field in the Yukawa model produces a scalar theory 
which has an infinite number of nonlocal interaction terms. Already taking only into account 
the first two orders of the 1/y expansion of the fermion determinant (see for example [|^) 
leads to a large number of non-single-site terms and also to terms which are no longer 
bilinear in the scalar fields. A small y expansion leads to infinite range interaction terms. 
The Yukawa model is therefore much more complicated than the metamagnetic Ising model 
and so its scalar propagators are analysed in this paper in a different way. In particular, a 
gaussian model does not describe the data. But the qualitative feature of the propagators, 
namely the existence of a second pole at ap = {it, it, it, vr) if two scaling regions overlap, can 
be understood considering simple gaussian and Ising models. 
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Figure Captions 



Fig. 1 Phase diagram of the 811(2)^® SU(2)r Yukawa model with naive fermions [|T^]. 
Fig. 2 Schematic description of the spectrum of scalar states in the phase diagram around 
point A. The staggered quantities are denoted by the index "st". We set a = l in the figure. 
Fig. 3 Schematic possible evolution of the bare Yukawa coupling constant i/b (its definition 
is described in the text) along a line of constant physics ?/i^=const in the FM(W) phase. For 
comparison also yo is shown. 

Fig. 4 Typical forms of the scalar propagators in the FM phase, mostly at negative k. 
Plotted is the inverse vr propagator in the momentum space as a function of ap^ = 2 Y^^ii^ ~ 
cos(ap^)). The left column shows results near the FM(W)-PMW phase transition and the 
right column near the FM(S)-PMS transition. With k decreasing downwards the points A 
and B, respectively, are approached. In this figure and in the following, error bars are not 
shown whenever they are smaller than the symbols. 

Fig. 5 An example of the inverse vr propagator in the FM(S) phase close to the point B. 
The full line is the fit with the two pole Ansatz analogue to (|6.3|) . 

Fig. 6 An example of the inverse vr propagator in the FM(W) phase at small positive k. We 
show in the upper part the Monte Carlo data and in the lower part the fit with the Ansatz 



Fig. 7 An example of the inverse vr propagator in the FM(W) phase at negative k close to 

the point A. We show in the upper part the Monte Carlo data and in the lower part the fit 

with the Ansatz ( |7.8| ) where a staggered scalar pole is also taken into account. 

Fig. 8 Examples of lattice size dependence of amp and avR. The dotted and dashed lines 

are fits linear in L~^. Extrapolations to infinite volume are also shown. 

Fig. 9 The wave function renormalization constant close to the FM(W)-PMW phase 

transition as a function of decreasing k. The point A is situated at k = —0.75(3). The 

dashed line serves to guide the eye. 

Fig. 10 Our present results for yn, as defined by the relation (|2.13| ), obtained from the 



infinite volume extrapolated values of amp and av^. The range of the earlier results |23 
obtained with a naive analysis of the Goldstone propagators and on finite lattices in indicated 
by the dotted error bars around the dotted curve (p.l|). 

Fig. 11 Ratio rria/vR as a function of the scalar correlation length The full squares 
represent a tentative infinite volume extrapolation. The circles are Monte Carlo results in 



the pure 0(4) $ theory |31 



Fig. 12 Phase diagram of an Ising metamagnet. The numerically determined positions of 
critical lines are indicated by the points with error bars and the dotted lines are mean field 
results. 

Fig. 13 Inverse momentum space propagator in the Ising metamagnet and its fit by means 
of the gaussian model with a nnn-term. 



